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Abstract 

We study the large N reduced model of D-dimensional Yang- Mills theory with special 
attention to dynamical aspects related to the eigenvalues of the N x N matrices, which 
correspond to the space-time coordinates in the IIB matrix model. We first put an 
upper bound on the extent of space time by perturbative arguments. We perform a 
Monte Carlo simulation and show that the upper bound is actually saturated. The 
relation of our result to the SSB of the U(1) D symmetry in the Eguchi-Kawai model 
is clarified. We define a quantity which represents the uncertainty of the space-time 
coordinates and show that it is of the same order as the extent of space time, which 
means that a classical space-time picture is maximally broken. We develop a I/O 
expansion, which enables us to calculate correlation functions of the model analytically. 
The absence of an SSB of the Lorentz invariance is shown by the Monte Carlo simulation 
as well as by the 1/D expansion. 
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1 Introduction 



Large N reduced models, which are the zero-volume limit of gauge theories, have been 
developed for studying the large iV limit of gauge theories J], ^, |3|, |, [5[]. Recently the 
(partially) reduced models have been revived in the context of nonperturbative formulations 
of string theory || The IIB matrix model ||, [|, which is the reduced model of 
ten-dimensional supersymmetric Yang-Mills theory, is expected to be a nonperturbative 
formulation of superstring theory. It is hoped that all the fundamental issues including the 
origin of the space-time dimensionality can be understood by studying the dynamics of this 
model. 

Aiming at such an end, we study the large N reduced model of Yang-Mills theory, which 
we call as the "bosonic model", since it is nothing but the bosonic part of the IIB matrix 
model. The action of the bosonic model is given by 

5 = -^tr([^,^] 2 ), (1.1) 

where (// = 1, • • • , D) are N x N traceless hermitian matrices. The partition function is 
given by 

Z = ( dA e" 5 , (1.2) 



where the measure dA is defined by 



D 



d.4 = n 



N / N 



;i.3) 



/i=l L*=l \i=l / i<j 

In Ref. [|1L)||, the value of the partition function of this model has been numerically 



evaluated and found to be finitef] for N > 4 when D — 3, for N > 3 when D = 4, and for 
N > 2 when D > 5, which means that in the large iV limit, the model is well-defined for 
D > 3 without any regularizations. This property, together with the fact that the action is 
homogeneous with respect to A^, enables us to absorb the coupling constant g by rescaling 
A'^ = -^A^. Hence g is nothing but a scale parameter like the lattice spacing in lattice 
gauge theories, and the dependence of expectation values on g is determined on dimensional 
grounds, which is also the case for the IIB matrix model ||. This is in striking contrast to 
the ordinary Yang-Mills theory before being reduced, where the parameter g is the coupling 
constant, which is independent of the scale parameter. 

: We will give an intuitive understanding for their conclusion in Section |2.1.l[ 
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We investigate dynamical aspects of the model related to the eigenvalues of A^, which 
correspond to the space-time coordinates in the IIB matrix model. One of the most im- 
portant quantities is the extent of space time defined by R = J(jftr(A 2 )). Since R should 
be proportional to yfg on dimensional grounds, we parameterize its large N behavior as 
R ~ ^fgN w . In the IIB matrix model, the value of to plays an important role when we 
deduce the space-time dimension, which is to be determined dynamically within the model. 
We determine the value of u for the bosonic model. We first obtain an upper bound on uj, 
namely uj < \ by perturbative arguments. The first argument is based on a perturbative 
calculation of the effective action for the eigenvalues of A^. The second one is based on exact 
relations among correlation functions derived as Schwinger-Dyson equations. We evaluate 
the correlation functions as a perturbative expansion and estimate the large N behavior of 
each term, from which we extract an upper bound on the extent of space time by requir- 
ing that the Schwinger-Dyson equations should be satisfied if we take all the terms of the 
perturbative expansion into account. We then perform a Monte Carlo simulation and find 
that uj = 7, which means that the upper bound given through perturbative arguments is 
actually saturated. This means that the perturbative estimation of the large N behavior of 
correlation functions is valid, which is also confirmed directly by the Monte Carlo simulation. 
We further clarify the relation of our result concerning the extent of space time to the SSB 
of the U(1) D symmetry in the Eguchi-Kawai model. 

The space time which is generated dynamically in the IIB matrix model is not classical, 
since A^s which are given dynamically are non- commutative generically. To what extent 
a classical space-time picture is broken is therefore an important issue to address. For this 
purpose, we define a quantity which represents the uncertainty of the space-time coordinates. 
We determine its large N behavior for the bosonic model and show that it is of the same order 
as the extent of space time, which means that the classical space-time picture is maximally 
broken in the bosonic model. 

In the IIB matrix model, the (ten-dimensional) Lorentz invariancef] should be broken 

spontaneously if the space time generated dynamically is to be four-dimensional. We define 

an order parameter for the SSB of the Lorentz invariance. We show, for the bosonic model, 

that the order parameter vanishes in the large N limit, and hence the Lorentz invariance is 

not spontaneously broken. 

2 When we define the IIB matrix model nonperturbatively, we need to make the Wick rotation and define 
the model with Euclidean signature. This is also the case for the bosonic model. Therefore, by Lorentz 
invariance, we actually mean the rotational invariance. 



2 



As a method complementary to the Monte Carlo simulation, we develop a 1/D expansion, 
which enables us to calculate correlation functions of the model analytically. To all orders of 
the 1/D expansion, we determine the large N behavior of correlation functions, which agrees 
with the one obtained by the perturbation theory and the Monte Carlo simulation. We also 
prove that the large N factorization property holds for the correlation functions. In order 
to confirm the validity of the 1/D expansion for studying the large N limit of the model 
for finite D, we perform explicit calculations of correlation functions up to the next-leading 
terms in 1/D and compare the results with the Monte Carlo data. We conclude that there 
is no phase transition at some finite D and that the 1/D expansion is valid for any D > 3. 

This paper is organized as follows. In Section |2], we focus on the extent of space time. 
We first put an upper bound on to, namely u < \ by perturbative arguments. We then 
show that the upper bound is actually saturated by a Monte Carlo simulation. We further 
clarify the relation of our result to the SSB of the U(1) D symmetry in the Eguchi-Kawai 
model. In Section [|, we show that a classical space-time picture is maximally broken in the 
bosonic model. In Section |], we develop a 1/D expansion, with which we determine the 
large N behavior of correlation functions. In Section |5], the absence of an SSB of the Lorentz 
invariance is shown by the Monte Carlo simulation as well as by the 1/D expansion. Section 
|6] is devoted to summary and discussion. 



2 The extent of space time 

2.1 Perturbative arguments 

2.1.1 Upper bound from the effective action 

We first start by decomposing An into its eigenvalues Aj„ and the angular part Vn as 



A, = V^V}, 



(2.1) 



where = diag(A M i, • • • , X^n) and is a unitary matrix. Since A^ is traceless, we have 
Si -Vi = 0- The effective action W(\) for the eigenvalues can be defined by 



-W(A) 



dV 



n ri(v - A*) 2 - w2) 



A 4 k i>j 



(2.2) 
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The VEV of an operator O can be written as 

JdXO(X) e" ff W 



(O) 



JdX e-^W 



(2.3) 



where 0(A) is defined by 
0(A) 



JdAjdV 


n,{rwv-^) 2 } 


6(A, - VpAtfj 


Oe~ s 


JdAjdV 


lit, {rii>j(v - Vj) 2 


} 8{A„ - V^Vi) 





(2.4) 



We can perform a perturbative expansion with respect to g 2 regarding X^i as external 
fields. The perturbation theory is valid when Aj's are widely separated, as we will see shortly. 

We first integrate out perturbatively around the unit matrix 3 either in ( |2.2| ) or in ( |2.4| ). 
We find that the van der Monde determinants are cancelled out and obtain the following 
constraints [01. 



d 



Q'fj.ijQ'liji 



0(a 3 ), 



(2.5) 



where d^i and a^ij are the diagonal parts and the off-diagonal parts of A^ respectively. 
Since the quadratic term in a of the action is 

1 



5, 



2g'- 



(2.6) 



we have to fix the gauge. Here we put the following gauge fixing term and the corresponding 
Faddeev-Popov ghost term. 

1 



S gJ- 

Sgh 



2g' 



? tr([A M ,a MJ 



-— tr([A M ,6][d M + a M ,c]) 



(2.7) 
(2.8) 



The total action S' = S + S g .f. + S g h can be written as follows. 



S' 
S' 



S 2 + S{ n t, 

y i-i-A v / 



(2.9) 
(2.10) 



Si 



hit 



tr([A pi a v ] [a M , a v \ + - [a M , a u ] 2 — b[X 



3 At first sight, one might think that should be expanded not only around the unit matrix but also 
around £W whose matrix elements are given as (E^)ij = S^r^j, where a is an element of iV-th order 



symmetric group. However, one can easily see that this does not affect the final results for (2.2) or (|2.4j 
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+ E {e 



„ A^j A^j 

7 Q'liikQ'fikiQ'vijQ'vji 



7 7 ® J [iikQ'itki\Q'iiij@'vji ~r Q>[iji"'vij ) 

+ V] ^ — T^a^ka^kiibijCji + bjiCij) \ + 0{a b 



(2.11) 



From ( |2.10| ), we can read off the propagators as 

((a^ija u ki)) = g 2 S^SuSjk 



(A, - A^ 
1 



(Cij6fci» = £ Watt rro, (2.12) 



where the symbol ((•••)) is defined by 

_ / dad&dc e~^ ^ 
f dad&dc e~ 5 2 

The interaction vertices can be read off from Si nt . 

The one-loop effective action W\{\) for A^ is given by 

Wi(A) = p - 2) J2 - ^) 2 ; ( 2 - 14 ) 

which is an 0(N 2 ) quantity and gives a logarithmic attractive potential^ between all the 
pairs of Aj. If the Aj's are sufficiently separated from one another, the one- loop effective 
action ( p. 14 ) dominates and higher loop corrections can be neglected. 

We can therefore examine the convergence of the integration over in (|1.2j ) in the 
infrared region, where the eigenvalues are far apart, by using the one-loop effective action 
(|2.14j) . The power of A coming from the one-loop effective action is — (D — 2)N(N — 1), 
while the one coming from the measure is D(N — 1). If and only if the sum of the powers is 
strictly negative, the integration over converges in the infrared region. Thus, we obtain 
the condition for the infrared convergence as 

iV>— — . (2.15) 
D-2 K ' 

Let us next discuss the behavior of the model in the ultraviolet region, where some pairs 

of the eigenvalues are close to each other. Looking at the one-loop effective action fl2.14D , 



As we review in Section 2.3, the bosonic model is actually equivalent to the weak coupling limit of 
the Eguchi-Kawai model. In this correspondence, the logarithmic attractive potential obtained here for the 
former model is essentially the one obtained 0] for the latter model. 
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one might argue that there might be an ultraviolet divergence. The two-loop effective ac- 



tion ( 2.23 ), which we calculate later, actually shows even a severer divergence for coinciding 
eigenvalues, and the higher loop one considers, the severer divergence one encounters. How- 
ever, this is simply due to the fact that the perturbative expansion is no more valid in the 
ultraviolet region, and does not imply a real divergence of the model. Rather, one should 
look at the partition function before performing the perturbative expansion: 



dA / d^ 



n{n(v-v) 2 } 



(2.16) 



where one sees no source of divergence for coinciding eigenvalues. 

Therefore the model is well defined if and only if the condition ( 2.15|) is satisfied. This is 
in agreement with the results obtained by a numerical evaluation of the value of the partition 



function |T0[. As far as the large N limit is concerned, the bosonic model is well defined for 
any D > 3. 







(a) 



(b) 



(c) 



(d) 



Figure 1: The diagrams we have to consider to obtain the two- loop corrections to the effective 
action for A M «. The wavy line corresponds to a^j and the dashed line corresponds to the 
6,c-ghosts. 



As we have seen above, if Aj's are widely separated, the Aj's attract one another due to 
the logarithmic attractive potential induced by the one-loop calculation, and the space time 
shrinks until the dominance of the one-loop effective action no more holds. This means that 
we can put an upper bound on the extent of the distribution of Aj by considering higher loop 
corrections to the effective action. 

Let us calculate the two-loop corrections. The diagrams we have to evaluate are shown 
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in Fig. [TJ The diagram (a) is evaluated as follows. 

(a) = ^M^a^^-a Ma a^ 

" \\ ijkl 



Defining 



4 2 



^■ui ^i/j ^ui ^uj 

Apt A/ifc <*\iii '•/ife 

& fj,ij & fij kj) ({&v klO'vlif)) 



ijkl 



+ y"UZ l^T^ r-^-((( a fMikO'nki)) {{ a vij a vji)) + ^fliiikO'vjif) ^ifi ijki&vijt)) 

— T r-^(((0'Liik a nki)) (( a ^ij a uji)) + ((O/xife&i/ji)) ((dfikidfiij))) 

5 (((O/iifeO/ifci)) (( a l ijiQ'vij)) + {(0'fiik a fiji)) ((Q'[iki a isij)}) 



Jl = JJ^-a^-a^ (2 ' 18) 



J 



2 



EnT^TlI' ( 2 - 19 ) 



13 ? JL* v - C (a, - - x k y ' (2 - 20) 



we can write the above result as 



(a) = g 2 \~\d{D - l)h + |l>(£> - l)/ 2 + 2(D - 1)J 3 | • (2.21) 
The diagrams (b), (c) and (d) can be calculated in the same way. The results are 

(b) = \g 2 {D-l)h, 

(c) = -\g 2 h, 

(d) = -2g 2 (DI 2 + I 3 ). (2.22) 



7 



Summing up all these results, we obtain the two-loop effective action as 

W 2 (X) = -((a) + (b) + (c) + (d)) = g 2 U(D - 2) 2 h - ±D(3D - 7)/ 2 - 2(D - 2)/ 3 

(2.23) 

In order to estimate the order of magnitude of the two-loop effective action, we note that 

#=^£*?V (2-24) 

Assuming that A,'s distribute uniformly in a D-dimensional ball with the radius R as in Ref. 
0, we can estimate the order of magnitude of I\ and I 2 as follows. 

r3 1 f aD„ id, a d ! - f N3 



r2 1 1 



Similarly, we can estimate the order of magnitude of I3 as O J . Therefore the order of 
magnitude of the two-loop corrections is O (^-^-^. By simple power counting, we can see 

that the n-loop corrections are of the order of O I A^ 2 ( ) J , in general. 

Since the one-loop result for the effective potential can be considered as an 0(N 2 ) quan- 
tity, the higher loop corrections can be neglected as long as R > ^fgN 1 / 4 . Thus, we can put 
an upper bound on R as R < y/giV 1 / 4 . 

2.1.2 Perturbative estimation of correlation functions 

In this section, we explain how to estimate the order in N of correlation functions of A 
through perturbative calculations. We will see that the perturbative estimation is valid if 
the upper bound on R given in Section |2.1.1| is actually saturated. 
Here we perform the integration over a M in ( |2.4j ) perturbatively as 

fdadbdcOe~ s ' 
°< A > = VdMce-*- • (2 - 27) 

An example of the operator O we consider is ti{A 2 ). We calculate O(X) as a perturbative 
expansion with respect to g 2 , which gives a loop-wise expansion. We then perform the 
integration over A for the result of each order of the perturbative expansion of O(X) with 
the Boltzmann factor e~ w ^ defined nonperturbatively by 
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As an example, let us consider the perturbative estimation of (-^tr((A 2 ) 2 )). 



The order of magnitude of the first and second terms are estimated as 0(R 4 ) and (g 2 N), 
respectively. We can easily see that the n-loop corrections are of the order of O ( R 4 ' 



■ i ■ v •» - y 1 j 

:> 1 1 .,,.,1 n i r2 J 

n\ 

in general. Therefore, if the upper bound on R obtained in Section |2.1.1| is actually saturated, 
i.e., R ~ ^fgN 1 / 4 , all the orders of the perturbative expansion give the same order in N and 
we obtain (^tr((A 2 ) 2 )) ~ g 2 N. 

We can generalize the above argument to the correlation functions written as 

itr(^^)ltr(A^4) • • • ^tx(A_^A . (2.29) 

Pi P2 PQ ' 

If the upper bound on R is actually saturated, the above quantity can be estimated as 

(^NVy, where P = E?=i^- 

In Section |2.2|, we calculate R through a Monte Carlo simulation and see that the upper 



bound on R is saturated and that the perturbative estimation of the large N behavior of 
correlation functions given above is indeed valid. 

2.1.3 Upper bound from Schwinger-Dyson equations 

In this section, we provide another way to extract information on the extent of space time R. 
The strategy is to derive exact relations among correlation functions of A through Schwinger- 
Dyson equations and to evaluate each of the correlation functions by the perturbative method 
described in Section |2.1.2| . By estimating the large N behavior of the correlation functions, 



we will see that the exact relations can be satisfied only if R < ^fgN 1 / 4 holds. 
We expand as 



N 2 -l 



A » = E K ^ ( 2 - 3 °) 



a=l 



where t a are the generators of SU(iV). We consider the following Schwinger-Dyson equations. 

= f dA^- a (tr(tVge- s ) , (2.31) 

= J dA^ (tr(t a A u A x A p )e- s ) . (2.32) 



These equations lead to the following exact relations, respectively: 



(i) -(tr([A„A v ] 2 )) = D(N 2 -l)g 2 , (2.33) 

(ii) - 9 2 ^ {( N ~Jf) (Map + <WM - ^<W<w} M^)> 

= (tr(A u A\A p [A T , A r ]])>. (2.34) 

We estimate the large N behavior of each of the correlation functions which appear in the 
above equations. 

Let us first consider eq. (i). We calculate the leading term of the l.h.s. as follows. 

- (tr[A M , A u f) l0 ~ ng -2{(tr[A M ,a„][A M ,aJ) - (tr[A M , a u ][\ u , a^])} 



2 S < ^^("V* _ -\uj)('Vi ~ ^^ij)((( l uji a uij 



= 2(D- l)N(N-l)g 2 . (2.35) 

The discrepancy between Q2.35Q and the r.h.s. of (i) is an 0(g 2 N 2 ) quantity. This should 
be compensated by higher loop corrections to the l.h.s. of (i). The typical next-leading 
contribution of the l.h.s. of (i) is obtained as follows. 

(2 ' 36) 

which is an 0( 9 ^ ) quantity. In general, the n-th subleading contribution is 0(g 2 N 2 {^M-) n ). 
If R > y/giV 1 / 4 , all the subleading terms are suppressed compared with the leading term and 
the discrepancy at the leading order cannot be compensated. We therefore obtain an upper 
bound R < y^iV 1 / 4 , which is the same as the one obtained in Section [2.1.1| . 

Let us next consider eq. (ii). The l.h.s. of (ii) can be evaluated to the leading order of 
the perturbative expansion as follows. 



i i r /..\ fading 

l.h.s. Of 11 ~ 



- g2 h { { N - w) {5 " Jxp + 5 " p5ux) " w 6 ^} A 7 ' (2 ' 37) 



The r.h.s. of (ii) can be evaluated as follows. 

(ti(A u A x A p [A T) [A^A T }})) 



10 



l0 ~ ng (tr(a v \x\ p [\ T , [a M , A T ]])) + (tr(A^a A A p [A T , [a M , A T ]] 
+ (tr(A J ,A A a p [A r , [a M , A T ]])) + (tr(A^A A A p [a r , [a p , A T ] 
+ (tT(a v \x\ p [\ T , [A M , a T ]])) + (tr(A^a A A p [A r , [A p ,a r ] 
+ (tr(A ly A A a p [A T , [A p ,a T ]])) + (tr(A^A A A p [a T , [A p ,a T ] 

= -/i {(iV - + M^a) - <WM ( E A * 

- 2(D - 2)g 2 ( ^ (V-^)^a^ \ (2 3g) 

The discrepancy between ( 2.37 ) and (|2.38|) which is dominant in the large N limit comes 



from the second term in ( 2.38|) , which is estimated as 0(g N R ). Sub-leading contributions 

9 2 N 
R 4 



to both sides of (ii) can be estimated as 0{g 2 N 2 R 2 {^-) n ). Thus we obtain the same upper 



bound on R. 

2.1 A Subtlety for the D = 3 case 

When we estimated the order in N of the expressions such as (|2.18|) , (|2.28| ) and (|2.36|) , we 
considered only the infrared region, in which the Xi S £1X6 widely separated from one another. 
Actually, we have to be careful about the contributions from the ultraviolet region, in which 
the Aj's are close to one another. As we will see below, the ultraviolet contributions do not 
affect the upper bound on R for D > 4, but they do for D = 3. 

Let us denote the minimum distance of the Aj's by r. The ultraviolet contributions can be 
estimated for ( p.!8| ), the second term of (|2.28| ), and (|2.36|) as ^, and s -^-, respectively. 



r 

.2 , 



For higher loop corrections, the above expressions should be multiplied by ) n , which means 
that the perturbative arguments break down for r < y/g. On the other hand, if r > ^/g, the 
ultraviolet contributions can be safely neglected in the perturbative arguments given in the 
previous sections. One should also consider mixed contributions such as the one obtained 



from ( 2 . 1 8|) when A, and Xj are near but Aj and A& are far apart, which can be estimated as 
Such contributions can also be neglected if r > ^fg. 
In order to estimate the magnitude of r, we introduce the typical distance between two 
nearest-neighboring Aj's, which is given byQ i ~ R- N^ 1 ^ . In general, we have r < £. If 



> 



y/g, the distance of the two nearest Aj's is controlled by an SU(2) matrix model 



5 If the Lorentz invariance were spontaneously broken, the D in the definition of I should be replaced by 
the dimension of the space time dynamically generated. We will see in Section [| that the Lorentz invariance 
is not spontaneously broken in the bosonic model. 
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and therefore, r > ^fg. Therefore, in order that the perturbative arguments may be valid, 
£ > ^/g, or equivalently, R > ^JgN l / D must be satisfied. Taking this point into account, the 
upper bound R < y^gN 1 ^ remains unchanged for D > 4, but for D = 3 case, it should be 
weakened to R < ^fgN 1 ^. 

2.2 Monte Carlo calculation of the extent of space time and other 
correlation functions 

In this section, we perform a Monte Carlo simulation of the bosonic model in order to 
determine the large N behavior of correlation functions. Above all, we determine the large 
N behavior of the extent of space time, defined by R = J {j^ti{A 2 )). The details of the 
algorithm for our simulation are given in the Appendix. 
The observables we measure are the following. 

(1) (^tr(A 2 )^x(^iV 1/4 )- 2 , 

(2) (^tr((^A,) 2 )^ x {VgN 1/A Y\ 



(3) ^-tr((A 2 ) 2 )^ x (VgN 1 ^) 

(4) (]^ 2 ))x (v^/y 4 , (2.39) 

where we define F^ v = i[A^, A v ]. Actually (4) can be written as (4) = 2 x ((3) — (2)). Note 
that the above quantities are independent of g, since we have normalized them so that they 
are dimensionless. 

We first note that <^tr(F 2 )) can be obtained analytically by the following scaling argu- 



ment. Rescaling as A^ — > e^A^ in the definition of Z in (L2), we obtain 
Z = e^ 2 - 1 ) f dA e"^ e4Etr(F2) 



= Z + ei ] D{N 2 - 1)Z- J dA tr (^ 2 ) e" 5 | +0(e 2 ). (2.40) 

Thus we obtain 

(^r(F 2 )^M(l-i) 3 2 (2.41) 



This agrees with the result ( |2.33| ) obtained through the Schwinger-Dyson equation. 

In Fig. |2], we show our results of the Monte Carlo simulation for D = 4 with N = 4 ~ 256. 
One can see that the quantities (1)~(4) are constant for N > 16. In order to see the finite 
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10' 
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0(1) 

□ (2) 
0(3) 

A (4) 



-B- ^ -6- 
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-3- 



-e- -e- 



10 



100 



1000 



Figure 2: The quantities (l)(circles),(2)(squares),(3)(diamonds) and (4) (triangles) are plot- 
ted against N for the bosonic model with D = 4. The solid line represents the exact result 
( |2.41| ) for (4). The tilted triangles are the plot of (A 2 ) /(y/gN 1 ^) 2 , which is discussed in 
Section [| 
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N effects, we plot in Fig. |3] the same data against 1/N 2 . The data can be fitted nicely 
to a Q + jik + jfc. In Section [4~2"1 , we will see that this large N behavior of the observables 
(1)~(4) agrees with the one given by the 1/D expansion. 



10.0 
9.0 
8.0 
7.0 
6.0 
5.0 
4.0- 
3.0 



2.0" 
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OO) 
A (4) 



ae 




( b coo q -----e- --- 



0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

1/N 2 



Figure 3: The quantities (l)(circles),(2)(squares),(3)(diamonds) and (4) (triangles) are plot- 
ted against 1/N 2 for the bosonic model with D = 4. The solid line represents the exact 
result (|2.41|) for (4). The other lines represent the fits to a + fk + 



In Fig. |], we plot the data for D = 10 with N = 2 ~ 32. Here again we find that the 
quantities (1)~(4) are constant for iV > 16. We have also checked that the data can be 
fitted nicely to a + p + p. 

It is rather surprising that the leading large N behavior shows up at about the same N for 
D = 4 and D = 10. If the system were a gas model with N particles in a D-dimensional space 
time, the large N behavior would be expected naively to show up at such iV that satisfies 
N 1 / ^> 1. In Section |^, we define the uncertainty of the space-time coordinates and show 
that it is of the same order as the extent of space time, which means that the bosonic model 
indeed can hardly be considered as a gas model with iV particles in a D-dimensional space 
time. 

From the results of the Monte Carlo simulation, we find that the extent of space time 
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Figure 4: The quantities (l)(circles),(2)(squares),(3)(diamonds) and (4) (triangles) are plot- 
ted against N for the bosonic model with D = 10. The solid line represents the exact result 
( [2.41Q for (4). The tilted triangles are the plot of (A 2 ) /(y/gN 1 ^) 2 , which is discussed in 
Section [| 
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R = ^(^tr(A 2 )) is of the order of O^iV 1 / 4 ), which means that the upper bound on R 
given in the previous sections for D > 4 is actually saturated. Due to this, the perturbative 
estimation of the large N behavior of correlation functions given in Section [2.1.2| with the 
input of R ~ 0(^/~gN 1 l A ) is valid for D > 4, which is also confirmed directly by the Monte 
Carlo simulation. 



10' 



10" 



• (1) D=3 
■ (2) D=3 

♦ (3) D=3 
A (4) D=3 

0(1) D=4 
□ (2) D=4 
O (3) D=4 
A (4) D=4 



■ ■ ■ 



▲ ▲ ▲ ▲ A ▲ 

^^eeeft # # # # 



10 



100 



1000 



Figure 5: The quantities (l)(circles),(2)(squares),(3)(diamonds) and (4) (triangles) are plot- 
ted against N for the bosonic model with D — 3. The results for the D = 4 are replotted 
for comparison. The solid line and the dashed line represent the exact result (|2.41|) for (4) 
for D = 3 and D — 4, respectively. 



We have done a simulation for the D = 3 case also, since it is rather an exceptional case 
from the perturbative point of view, as we remarked in Section |2.1.4| . The result is shown in 
Fig. ^ One can see that the result is qualitatively the same as for the D = 4 case. The fact 
that the upper bound R < ^J~gN 1 ^ is not saturated means that the perturbative arguments 
completely break down in the D = 3 case due to the ultraviolet contributions. We will see 
in Section [4.3| , however, that the 1/D expansion is valid even for D = 3, which explains why 
the results for D = 3 are qualitatively the same as those for D > 4. 
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2.3 Relation to the U(1) D SSB of the EK model 

In this section, we examine the relation of our result concerning the extent of space time to 
the U(l) £) SSB of the D-dimensional Eguchi-Kawai model. 

The action of the D-dimensional Eguchi-Kawai model is given by 

D 

S EK = -N/3 HU»U v Upt). (2.42) 
This model has a U(1) D symmetry. 

y M -e*"EV (2.43) 

In Ref. P], it has been shown that if the U(1) D symmetry is not spontaneously broken, the 
model is equivalent to an SU(iV) lattice gauge theory on an infinite lattice in the large N 
limit, where the coupling constant /3 in the action ( |2.42|) is kept fixed. 

It is found in Ref. that when D is larger than two, the U(1) D symmetry is sponta- 
neously broken in the weak coupling (large f3) region. This means that we cannot study the 
continuum limit of the large N gauge theory by the EK model when D > 2, since we have 
to send (3 to infinity when we take the continuum limit. Slightly modified versions of the 
model known as the quenched EK model [Q, f|, [| and the twisted EK model [T1J have been 
shown to be equivalent to the SU(iV) lattice gauge theory on the infinite lattice in the large 
N limit. 

Here we note that the bosonic model is actually equivalent to the weak coupling limit of 
the original EK model for D > 2. This can be seen as follows 0. Since is unitary, the 
eigenvalues of are on a unit circle in the complex plane. The U(1) D symmetry rotates all 
the eigenvalues by the same angle. In the weak coupling limit, namely when /3 — ► oo, the 
U(1) D symmetry is spontaneously broken and the eigenvalues collapse to a point. Therefore, 
when p is sufficiently large, the dominant configurations are given by 

C/ M ~ e^e^, (2.44) 

where a M are real numbers defined modulo 27r and are NxN traceless hermitian matrices, 
whose eigenvalues are small. a M take random values due to the U(1) D symmetry. The action 
for is obtained by putting ( [2.44j ) in ( 2.42 ) and expanding it in terms of A^. We obtain 
the following. 

1 D 

Sek ~ ~2 N P X] tr (tA"^M 2 ) + (higher order terms in A). (2.45) 
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In the weak coupling limit, the higher order terms can be neglected and we arrive at the 
bosonic model fll.lj) with g = ^L^ . 

A quantity which represents the extent of the eigenvalue distribution in the Eguchi-Kawai 
model can be defined by || 



As can be seen in the above expression, < P < | for any U^. P ~ | when the eigenvalues 
are uniformly distributed, while P ~ when the eigenvalues collapse. Thus (P) serves as 
an order parameter for the SSB of the U(1) D symmetry. 

Let us consider (P) in the weak coupling limit. Putting ( |2.44| ) in ( p.46| ) and expanding 
it in terms of A^, we obtain, to the leading order of A^, 

Therefore (P) ~ ^R 2 i where R is the extent of the eigenvalue distribution defined in the 
bosonic model. In Section 2^2, we have seen that R 2 = C{^fgN 1 ^) 2 in the large N limit. 
This means that 

<p > ~ ^BW' (2 - 48) 

as (3 —>■ oo. Thus our finding for the bosonic model can be interpreted in terms of the EK 
model as the eigenvalue distribution having a finite extent in the large N limit for fixed (3 
in the weak coupling region. We check this explicitly by a Monte Carlo simulation in what 
follows. 

A Monte Carlo simulation of the EK model with D = 4 and N = 20 has been performed 



in Ref. [fCql , where the internal energy defined by 

E = ^^^Etr(W^)) (2-49) 
is measured. The asymptotic behavior of the quantity in the weak coupling limit can be 



obtained just in the same way as we derived ( [2.48] ). The result to the leading order of A^ 
reads [112 1 



E ~ 1 + 1- -/ S^tr(\A u ,A 

2ND D - 1 \^ u M ' 
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Figure 6: The order parameter (P) for the U(l) r> symmetry (open squares) as well as the 
internal energy E (filled squares) is plotted against (3 for the D = 4 SU(16) Eguchi-Kawai 
model. The solid line and the dashed line represent the asymptotic behavior of (P) and E 
respectively at large (5 predicted by the bosonic model. 
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where we used the exact result (2.41) in the second equality. 

In Fig. H we plot the order parameter (P) as well as the internal energy E for the SU(16) 
EK model. The solid line and the dashed line represent the asymptotic behavior at large (3 
for (P) (p.48|) with C = 2.162(5) and that for E (|2.50|) , respectively, predicted by the bosonic 
model. The data for large (3 fit nicely to the predictions. The results for the internal energy 
are in good agreement with the results in Ref . |[TjJ , though N differs from ours slightly. We 
have also done simulations with N = 32 and N = 64 for (3 = 1.0 and checked that both E 
and (P) remain the same within error bars. This large N behavior is just the one we expect 
from the large N behavior of the quantities (1) and (4) in the bosonic model obtained in 
Section f2~2| . 

(P) is around 0.5 for (3 < 0.17 and decreases for (3 > 0.17, which shows that the U(1) D 
symmetry is spontaneously broken at around f3 = j3 c ~ 0.17. At the critical point (3 = /3 C , 
the internal energy is continuous but its first derivative with respect to (3 seems to diverge, 
which suggests that the system undergoes a second order phase transition. 

It is intriguing to compare the phase diagram of the EK model with that of the TEK 
model [11, HJ. Since the U(1) D symmetry is not spontaneously broken in the TEK model 
even in the weak coupling region, the model is equivalent to the SU(iV) lattice gauge theory 
on the infinite lattice in the large N limit throughout the whole region of (3 ]TI] . The system 
undergoes a first order phase transition at (3 ~ 0.35. This phase transition can be considered 
as an artifact of the Wilson's plaquette action, and could be avoided if one wishes, say, by 
introducing the adjoint term . 

If the (3 C , at which the SSB of the U(1) D symmetry occurs in the EK model, were larger 
than 0.35, the EK model would also have undergone a first order phase transition at (3 ~ 0.35, 
but the fact is that (3 C < 0.35 and the model is not equivalent to the SU(iV) lattice gauge 
theory for f3 > j3 c . Thus the first order phase transition may well be absent in the EK model. 



Before ending this section, let us address the issue [10] whether the bosonic model can be 
used to calculate correlation functions in the large N Yang-Mills theory. It is natural to think 
that the answer is in the negative considering the fact that the EK model is not equivalent 
to the large iV lattice gauge theory for D > 2 due to the SSB of the U(1) D symmetry. Let 
us see how this can be restated within the bosonic model. 

One way to show the equivalence between the reduced model and the original model in 
the large N limit is to confirm it to all orders of the perturbative expansion [^, [J, where 
X^i should be identified with the loop momenta and the extent of space time R plays the 
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role of the momentum cutoff. One can easily see that the coupling constant gyu of the 
corresponding Yang-Mills theory is given by 

9ym = ^p, (2-51) 

which can be also deduced on dimensional grounds. In order to take a nontrivial large 
iV limit, namely the 't Hooft limit, one has to fix g^M-^ = Xym to a constant. Since 
R = C 1 ! 2 ^fgN 1 ^ , where C is a fixed constant determined dynamically, we obtain 

\ YU = C-%(g 2 N) 1 -T. (2.52) 

From fl2.52j ), we have to fix g 2 N to reproduce the 't Hooft limit. In this limit, R goes to a 
constant, which we denote by A. Xym can then be written as 

Xym = C~ 2 A 4 ~ D , (2.53) 

which means that we are not allowed to take the continuum limit A — > oo as one wishes, but 
the coupling constant Xym is doomed to scale canonically with the cutoff A. Thus, there is 
no way to take a nontrivial continuum limit in the Yang-Mills theory to which the bosonic 
model is equivalent perturbativelyQ. We therefore conclude that the bosonic model cannot 
be used to calculate correlation functions in the large N Yang-Mills theory. 



3 Breakdown of the classical space-time picture 

As we mentioned in the Introduction, the eigenvalues of represent the space-time coordi- 
nates, when we interpret the IIB matrix model as a string theory formulated nonperturba- 
tively. If A^s are mutually commutative, they can be diagonalized simultaneously and the 
diagonal elements after diagonalization can be regarded as the classical space-time coordi- 
nates. This must be violated more or less by A^s generated dynamically in the IIB matrix 
model. Therefore it makes sense to ask to what extent the classical space-time picture is 
broken. In order to address this issue, we define a quantity which represents the uncertainty 
of the space-time coordinates. We determine its large N behavior for the bosonic model and 
show that it is of the same order as the extent of space time, which means that the classical 
space-time picture is maximally broken in the bosonic model. 

6 The equivalence does not hold in the strict sense even perturbatively. See Ref. 0] for the details. 
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We define such a quantity by considering an analogy to the quantum mechanics. We 
regard the matrices A^s as linear operators which act on a linear space, which we identify 
as the space of states of particles. We take an orthonormal basis |ej) (i = 1, 2, ■ • • , N) of the 
N- dimensional linear space, and identify the ket |e$) with the state of the z-th particle. The 
space-time coordinate of the i-th particle can be defined by {e^A^i). The uncertainty 5(i) 
of the space-time coordinate of the i-th particle can be defined by 

^) 2 = E{^i^ 2 i^)-(^i^i^)) 2 }- (3- 1 ) 

Note that this quantity is invariant under a Lorentz transformation and a translation : 
A^ — > A^ + a M . We take an average of S(i) 2 over all the N particles. 

i 

Note that this quantity depends on the orthonormal basis |e$) (i = 1, 2, • • • , N) we choose. We 
therefore define a quantity A which represents the uncertainty of the space-time coordinates^ 
for a given A^ by 

A 2 = min{¥) (3.3) 

{>;}} 

= ^ti{A 2 )- max ±-Y{(UAM%} 2 . (3.4) 

i 

Note that A is invariant under a gauge transformation: A^ — > gA^g*. Note also that A = 
if and only if A M are diagonalizable simultaneously. The classical space-time picture is good 
if A is smaller than the typical distance £ ~ R • iV _1//D between two nearest-neighboring 
particles. 

We calculate (A 2 ) for the bosonic model by a Monte Carlo simulation. Given a configu- 
ration A At , which is generated by a Monte Carlo simulation, we evaluate A 2 numerically by 
performing the maximization of ^2i{{U A^U ')j«} 2 in fl3.4p in the following way. We first max- 
imize ^2i{{U A^U^)n} 2 restricting the U to be in one of the SU(2) subgroup of the SU(iV). 
Redefining the A^ by the UqA^Uq, where Uq is the element of the SU(2) subgroup that gives 
the maximum, we do the maximization for all the SU(2) subgroups successively. We 

call this as 'one sweep' in our maximization procedure. For N < 32, the quantity saturates 

7 It might be tempting to identify the quantity A with the space-time uncertainty which appears in 
the space-time uncertainty principle proposed in Ref. [ p"4| |. However, it is unclear at present whether this 
identification is correct. 
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up to 16 digits within 200 sweeps for D = 4, and within 500 sweeps for D = 10. For D = 4, 
we have done the measurement with larger N as well. For iV = 256, we have made 300 
sweeps of the maximization, with which the saturation is achieved up to 3 digits. The result 
thus obtained is plotted in Figs. |2] and f|. We find that (A 2 ) normalized by (^fgN 1 / 4 ) 2 tends 
to a constant for large N. Thus we conclude that a/ (A 2 ) is of the same order as R. 

We have seen in Section 7 TJ1 that the N at which the leading large N behavior dominates 
does not depend much on D. The large N behavior would show up for N 1 ^ ^> 1 if the 
system can be viewed as N classical particles in a D-dimensional space time. Indeed our 
conclusion in this section means that the bosonic model can hardly be viewed as consisting 
of N particles. 



4 1/D expansion of the bosonic model 
4.1 The formalism 

In this section, we show that the bosonic model allows a systematic 1/D expansion, which 
serves as an analytical method complementary to the Monte Carlo simulation. To all orders 
of the 1/D expansion, we determine the large N behavior of correlation functions and prove 
the large N factorization property. We also perform an explicit calculation of the quantities 



(1)~(4) defined in Section ^]2| up to the next-leading order and compare the results with the 
Monte Carlo data for various D, from which we conclude that the 1/D expansion is valid 
down to D = 3. 

Here we use again the adjoint notation defined through 

N 2 ~l 
a=l 

where t a are the generators of SU(iV). Taking the generators t a to satisfy the orthonormal 
condition 

tr{tn b ) = 5 ab , (4.2) 

we have the following relation 

^2(t a )ij(t a ) k i = Su5 jk - —SijSki- (4.3) 
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S = ~— 2 tr([A„A u } 2 ) 



We first rewrite the action ( |1 . 1| ) for the bosonic model as 

1 

1 , 

where we have defined 

K abcd 1 Jx„/uo iCUj-b ■ ' ' a ^ b 



\ abcd A^A b A v c A d , (4.4) 



A abcd = S jtr([*V c ]M) + (a <- 6) + (c<-> d) + I ^ 1 j . (4.5) 

Note that the measure ( |1.3| ) can be written as JT ^A® up to an irrelevant constant factor. 
We introduce the auxiliary field h ab , which is a real symmetric tensor, with the following 
action. 

S[A,h] = -^ 2 \ abcd {h ab h cd - A*A b h cd - h ab A v c A d ) (4.6) 

= ^\^h ab h cd + ±K ab A;A b , (4.7) 

where we have defined the dimensionless kernel 

K ab = --\ abcd h cd (4.8) 
9 



for A a . Integrating out the auxiliary field, we reproduce the original action (O ). 

Since the A is quadratic in ( |4.7| ), we can integrate it out first. The propagator for A° 
is given by 

{A;aJ>) = gS^K-'U. (4.9) 
For example, (tr(A 2 )) can be expressed by the following integral. 



l2AN _ / &h ab gD(K~ l ) aa e- s ^ 



S e s is the effective action for h defined by 

S cS = ^TrlnK + -L\^ d h ab h cd , (4.11) 

where the Tr represents a trace over the adjoint indices a. 

By rescaling the h as h ab = —^h ab , which is now dimensionless, we can rewrite the 
effective action as 

S cS = j |Tr In K + h abcd h ab h cd \ . (4.12) 
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Therefore, in the large D limit, the integration over h is dominated by the saddle point and 
one can perform a systematic 1/D expansion. The saddle point equation can be given as 

(K~ l )J^ + X abcd h cd = 0. (4.13) 
oh ab 

The relevant saddle point h^ b , if we assume that it preserves the SU(iV) symmetry, should 
be written as = vd ab . From the saddle point equation ( |4.13j ), we obtain v = -j^- We 
expand the h ab around the saddle point h^} as 

i fW 

h ab = -^=5 ab + 2\ —ip abl (4.14) 



f2N V D 



where ip ab is real symmetric. We have put the factor 2^/yj in front of ip ab for later convenience. 
Expanding K in terms of <p, we obtain 



K ab = V2ND(5 ab - eX ab ), (4.15) 



where e = J ^ and 



X ab = X abcd <p cd . (4.16) 
We can expand (K~ l ) ab and Trln.fr in terms of ip in the following way. 

{K- l ) ab = -^==(5 ab + eX ab + e 2 X ac X cd + 0{^)), (4.17) 

TrlnK = const. - e X aa -^e 2 X ab X ba + 0{^). (4.18) 

The effective action S e s can be expanded in terms of ip as 

1 00 ,n-2 

S cS = -T abcd ip ab if cd - V Tr(X n ) + const., (4.19) 

2 ' n 

n=3 

where we have defined 

rpabcd yabef ycdef _|_ c^j^ yabcd 20) 

The propagator of ip ab can be given by 

Wed) = S abcd , (4.21) 

where S abcd is defined by 

rpabpq gpqcd gabpqrppqcd jabcd 22) 
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In order to obtain S abcd explicitly, we introduce the following quantities. 



jpabcd 


^ jabcd _|_ jbacd _|_ jabdc _|_ jbadc 


Qabcd 


^ ^acbd _|_ jadbc^ 


jjabcd 


VabOcdi 


jabcd 


l(5 ac 5 bd + 5 ad 5 bc ), 



where 

jabcd = ^a^dy 

\ abcd and T abcd can be written as follows. 



yabcd ^^Qaoca pa 

rpabcd _ QJ\J p abcd -|_ 4iVG a ^ C ^ <2Jjabcd ^jabcd 



yabcd T?abcd\ 



gabcd can chained as 



S' 



abed 



A(N 2 - 1] 



2 1 

rpabcd i ]\TS-iabcd , rjabed , jabcd 

3N F +NG + 6W H +I 



(4.23) 

(4.24) 
(4.25) 
(4.26) 

(4.27) 

(4.28) 
(4.29) 



(4.30) 



V 



abed 



Figure 7: 



Since we work with X ab instead of ip ab in what follows, we introduce the propagator of 
X ab given as 



V 



abed def 



(XabXcd) 

yobef gefgh^ghed 

1 f 2 



1 



j\t rpabcd , j\j s~iabcd , rrabed , jabcd 

N 2 - 1 1 3 6 



(4.31) 
(4.32) 
(4.33) 



Let us represent V abcd by the four-point vertex, as is shown in Fig. |7|. Then we can 
calculate (tr(A 2 )) by evaluating the diagrams depicted in Fig. §. The result is given as 
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o ,o, o 



(a) (b) (c) 

Figure 8: 



follows. 

(tr(A 2 )) = qB 1 { (N 2 - 1) + e 2 (V abab + V aabc V bddc ) + 0(e 4 ) } , (4.34) 
y2ND 

where each term corresponds to the diagrams (a), (b) and (c) in Fig. ^ respectively. By 
noting that V abcc = —\5 a b and V abab = |(7iV 2 — 1), we obtain the result for the quantity (1) 
defined in Section ET2| as 



;i) = (^tr(A 2 )^x(^iV^)- 2 

f{0-*H(™M*)}- ™ 

Similarly, we can calculate (tr^A^A^ 2 )) as follows. 

(tr((A,A u ) 2 )) 
= G abcd (A;A b A u c A u d ) 

= g 2 G abcd {D 2 ({K-% b {K-\ d )+2D{{K-% c {K- 1 ) bd )} 

g2 _ jjj^ r (jaabb | ^2 abed y abed _|_ 2^! aa ^ c ^^ c ' c _|_ ^Qabccyabdeydf fe\^ ~^ 

+ 2DG abab + 0{l)}. (4.36) 

(tr((y4 2 ) 2 )) can be calculated by replacing G by F in the above expressions. We obtain the 
following final results for the quantities (2) and (3). 

def , x 2 A v r/77 Ar l/4>-4 



(2) = ( jjtr((A,A u ) 2 ) ) x (^TgN^) 
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+ Ul-^ + &) +o &»> (4 - 37) 



(3) = 0((^) 2 )>x(v^^)- 



" +V\+o(-k)}. (US) 



D V2 2N 2 iV 4 J \D 2 
Note that the above results are consistent with the exact result (|2.41 ) for (^ti(F 2 )) . 

4.2 Large N behavior and factorization 

We estimate the order in N of correlation functions to all orders of the 1/D expansion and 
further prove the large N factorization property. 

Let us first consider (jjti(A 2 )) as an example. The diagrams we encounter up to the 
next-leading order in the 1/D expansion are shown in Fig. |j. We estimate the order in N of 
each diagram we encounter at each order of the 1/D expansion. We start by replacing the 
yabcd r.h.s. of ( [4.33 ). Each ^-vertex should be replaced by either of the F-, G-,H- 



and /-vertices with the factor for the first two, and with the factor -j^- for the last two 
due to the coefficients in (|4.33[ ). We denote the number of F- and G- vertices by V\ and the 
number of H- and /-vertices by V2 in the diagram after this replacement. The diagram has 
an overall factor (j^) Vl {^2)^ ' ■ 

The next task is to replace the F-,G-,H- and /-vertices by their definitions ([4.23 )^( 4.26 ) . 



Let us introduce the diagrammatic representation for tr(t a t b ■ ■ ■ t c ) as in Fig. |9[ Fig. [H] shows 



tr(t a t b -t c ) 




Figure 9: 



an example of what we get from the diagram (c) of Fig. |8|. The symbols F and H in the 
l.h.s. of Fig. [H] represent that the ^-vertices in the diagram (c) of Fig. ||have been replaced 
by F- and H- vertices respectively. In general, we obtain diagrams composed of V\ blobs 
connected by L lines with additional i loops which are not connected to any of the blobs. 
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H 





Figure 10: 



Reflecting the fact that each blob has four legs, we have 



AV l = 2L. 



(4.39) 



Let us denote by C the number of connected parts in the diagram excluding the additional 
£ loops. For the diagram shown in the r.h.s. of Fig. [It], we have C = 1 and £ = 1. Since 
the number of connected parts has a chance to increase from its initial value '1' by the use 
of H- and I- vertices, we have 

C + £ = 1 + j, (4.40) 

where j — 0, 1, • • • , V%. The £ loops give a factor of (N 2 ) e . 

In order to read off the order in N of the diagram, it is convenient to switch at this stage 
to the double-line notation. The tr(t a t b ■ ■ ■ t c ) should be represented as in Fig. O. Due to 



Figure 11: 



the relation ( |4.3j ), we have the diagrammatic rule depicted in Fig. [12] for the contraction 
of the adjoint indices. In order to see the leading large N behavior, it suffices to consider 
only the first term. Fig. ^ shows the diagram in the double-line notation which gives the 
leading large iV contribution to the diagram on the r.h.s. in Fig. |l(]. If we denote by F 
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J_ 

N 



Figure 12: 





Figure 13: 



the number of index loops excluding those coming from the I additional loops, we have the 
following relation 

F + V x - L = 2C - 2h, (4.41) 

where h is the number of handles necessary to write the diagram without crossings of the 
double lines. 

The order of magnitude of the contribution to (-^tr(v4 2 )) from the diagram is given by 

1 9 



lY ( J-)\ N zy N F 



(4.42) 



where we have used Q4.39| ),( fO0] ) and ( |4.41| ). Therefore the leading term is O(gyN), which 
is obtained for j = V2 and h — 0. 

So far, we have been concentrating on the leading large N behavior. Let us next consider 
sub-leading terms. Note first that the sub-leading terms from ( f4.42|) are given in terms of 
a 1/N 2 expansion. The effect of the second term in Fig. [12| is a factor either of 0(1) or of 
0(^2) for each replacement of the first term by the second term. Also the prefactor Af2 1 _ 1 
of ( |4.33| ) can be expanded in terms of 1/N 2 , which may give extra 0(jj%) factors, as well. 
Thus we have shown that the subleading terms appear as a 1/N 2 expansion. 
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Generalization of the above argument to the quantity written as 



PI P2 PQ 



(4.43) 



is straightforward. (|4.39|) , (|4.40|) and ( [4.4 1|) should be replaced byf] 



W 1 + P = 2L, 
C + £ = q + j, 
F+(V 1 + Q)-L = 2C-2h, 



(4.44) 
(4.45) 
(4.46) 



where P = Yli=i Pi anc ^ Q , = 1)2,---,Q. g is the number of connected parts before replacing 
the V^-vertices by either of F-,G-,H- and /-vertices. Then the order of magnitude of the 
quantity (|4.43| ) is given by 

P/2 / 1 \ Vi / \ V 2 



9 



IV' f J. 

( ^jN 1/4 ) p N 2{{c >~ Q)+{j ~ V2) ~ h \ 



(4.47) 



where we have used ( |4.44|) , (|4.45|) and ( |4.46| ). Therefore, we have the maximum order 
{^/gN l / 4 ) p when q = Q, j = V 2 and h = and the sub-leading terms appear as a 1/N 2 
expansion. 

Thus we have shown, to all orders of the 1/D expansion, that if we define 



(4.48) 



correlation functions of 0's such as (O1O2 ■ ■ ■ O n ) have finite large N limits. Note that this 
large N behavior is the same as the one obtained by the perturbative argument in Section 
2.1.2\ with the input of R ~ ^/gN 1 ^. The fact that the large iV behavior of the correlation 
functions does not depend on the order of the 1/D expansion is consistent with the fact 
that the perturbative argument is independent of D except for D = 3. We have also shown 
that the subleading large N behavior of the correlation functions can be given by a 1/N 2 
expansion, which is clearly seen by the Monte Carlo simulation in Section |2.2| . Moreover, 
the fact that the leading large N contribution to ( [4.43 ) comes solely from the diagrams with 



8 (fa\44|), J4.45D and Q4.46| ) do not reduce to fl4.39| ), (|4.40| ) and flOl] ), when we take Q = 1, px = 2. This is 
simply because we have omitted the blob coming from the trace of (tr(A 2 )) using the fact that ti(t a t b ) — 6 a b- 
The general argument below holds for the particular case as well. 
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(0,0, ■ ■ ■ O n ) = (O t )(0 2 ) ■ ■ ■ (O n ) + Ol±-) (4.49) 



q = Q as is seen from ( |4.47| ) means that we have the factorization property : 

TP) 

in the large N limit. 

The factorization can be generalized to the case when (9's are Wilson loop operators such 

as 

O = -^Vexp(i J dak^(a)A^). (4.50) 
We here consider k^(a) to be written as 

*» = -^irM^ ( 4 - 51 ) 

where /^(cx) is independent of g and N. Then the above quantity can be expanded in terms 
of / as 

oo „i -i „i 

= X/ " / d(J i / d(7 2 • • • / da n f^ («7i)/ w (02) • • • / Mn M 

n=0 ® 1 "Cn—l 

ItrOV • • -A M2 A m ) x (^jWt" (4.52) 

Therefore the factorization property is satisfied by the Wilson operators considered above 
as well. This is in contrast to the situation seen in the double scaling limit of the 2D 



Eguchi-Kawai model, studied as a toy model of the IIB matrix model |15 



4.3 Comparison of the 1/D expansion with numerical data 

In order to confirm that the 1/D expansion developed in the previous sections is valid for 
studying the large iV limit of the model for finite D, we perform Monte Carlo simulations 
for various D and examine the D dependence of the quantities (1)~(4) defined in Section 
|2.2| . Extrapolation to iV = 00 has been done by the fit to a® + ai/N 2 with iV = 16, 32, 
64, 128, 256 for D = 3, 4, with N = 16, 32, 64 for D = 6, 8, and with N = 16, 32 for 
D = 10,12,14,20. In order to compare the data with the results of the 1/D expansion 
( |4.35| ), Q4.37D and ( |4.38| ), where we put iV = 00, we normalize the Monte Carlo data by the 



leading asymptotic behavior in the large D limit. For (2), since the leading term in 1/D 
vanishes in the iV — >• 00 limit, we normalize the data by the next-leading term. In Fig. Ill 
we plot the results against 1/D. We fit the data to 1 + b t /D + b 2 /D 2 + b 3 /D 3 + b 4 /D 4 , 
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Figure 14: The quantities (l)(circles),(2)(squares),(3)(diamonds) and (4) (triangles) extrap- 
olated to N = oo and normalized by the leading asymptotic behavior in the large D limit 
are plotted against 1/D for the bosonic model. The solid line represents the exact result 
( ET41D for (4) and the other lines represent the fits to 1 + h/D + b 2 /D 2 + b 3 /D 3 + b A /D 4 with 
the input of the 1/D expansion up to the next-leading terms. 
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except for (4), which is known exactly. For (1) and (3), we use the values of b\ predicted by 
the next-leading term. We find that the data are nicely fitted. 

This result suggests that we do not have a phase transition at some intermediate value 
of D and that the 1/D expansion is valid for finite D down to D = 3. 



5 No SSB of Lorentz invariance 

In this section, we address the issue of the SSB of the Lorentz invariance in large N reduced 
models. This is of paramount importance in the IIB matrix model, since if the space time is 
to be four- dimensional at all, the 10D Lorentz invariance of the model must be spontaneously 
broken. 

Let us define an order parameter for the SSB of the Lorentz invariance. We first note that 
the extent of the eigenvalue distribution in the direction n^, where is a .D-dimensional 
unit vector, can be given by the square root of 

J(n) = -^tr(n-A) 2 (5.1) 

= ( 5 - 2 ) 

where is a D x D real symmetric matrix defined by 

V = itr^VU (5.3) 

The eigenvalues of are real positive. The SSB of the Lorentz invariance can be probed 
by the variation of the eigenvalues, which is given by 

2 



J = - (^I„ ] (5.4) 



If (J)/i? 4 is nonzero in the large N limit, where R ~ ^/gN 1 ^ is the (averaged) extent of 
space time, the Lorentz invariance is spontaneously broken. Therefore, ( J) / { y /gN 1 ^) i can 
be considered as an order parameter for the SSB of the Lorentz invariance. 



In Fig. [15] we plot ^^J^l /4 ^ against N for D = 4, 6, 8 and 10. We see that the order 



parameter vanishes in the large N limit, which means that the Lorentz invariance is not 
spontaneously broken. 
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Figure 15: The order parameter for the SSB of the Lorentz invariance (J) /(y/gN 1 ^ 4 ) 4 is 
plotted against iV for D = 4, 6, 8 and 10. The solid line represents the leading large N 
behavior for D = oo and the other lines represent the fits of the data to the 1/N 2 behavior 
predicted by (|4.49|) . 
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This result can be also confirmed by the 1/D expansion. By using the 1/D expansion, 
we can prove a relation 

HA,A V )) = l^(tr(A 2 )), (5.6) 

which already suggests the absence of an SSB of the Lorentz invariance. In fact, this relation 
combined with the relation 



Ltr(A„A„)±tr(A x A p )^ x (v^V^)" 4 



(5.7) 



due to the large N factorization property ( 4.4S ), which can be also shown by the 1/D 
expansion, leads to ^^\ /4 ^ = 0(jp) in the large N limit. Indeed the large N behavior of 
the order parameter obtained by the Monte Carlo simulation can be nicely fitted to I/TV 2 . 

We have also calculated the order parameter by the 1/D expansion explicitly. The result 
is given as 



(V^V 1/4 ) 4 2iV 2 [ 3D \D 2 J J V^ 4 
The leading large N behavior for D = oo is plotted in Fig. [15] for comparison. 



6 Summary and Discussion 

In this paper, we studied the reduced model of bosonic Yang-Mills theory with special 
attention to dynamical aspects related to the eigenvalues of the N x N matrices, which are 
regarded as the space-time coordinates in the IIB matrix model. We found that the bosonic 
model allows a systematic 1/D expansion, which indeed revealed many interesting dynamical 
properties of the model. One should note again that the parameter g of the model, which 
corresponds to the coupling constant before being reduced, is merely a scale parameter which 
can be scaled out if one wishes by an appropriate redefinition of the dynamical variables A p . 
Thus, unlike in ordinary gauge theories, we have neither a weak coupling expansion nor a 
strong coupling expansion. 1/D is the only parameter which allows a systematic expansion 
of the model. We confirmed that the 1/D expansion is indeed valid for all D > 3. 

To all orders of the expansion, we estimated the large iV behavior of correlation functions, 
which agrees with the one determined by the Monte Carlo simulation. Above all, the extent 



36 



of space time, which is defined by R = J {j^ti(A 2 )), was found to behave as R ~ ^fgN 1 ^ 
in the large N limit. For D > 4, the above result turned out to saturate the upper bound 
given by perturbative arguments, which fact makes the perturbative estimation of the large 
N behavior of correlation functions valid. We also showed that the large N behavior of the 
extent of space time is consistent with the SSB of the U{1) D symmetry in the weak coupling 
region of the Eguchi-Kawai model. 

The fact that we could obtain an upper bound on R through a perturbative evaluation 
of the effective action has an important implication. It means that if we look at the effective 
action for the eigenvalues obtained perturbatively, it has nonnegligible contributions 
when a pair of Aj's which are assigned to adjacent index loops on a Feynman diagram (e.g., 
Fig. [I]) are separated widely in the D- dimensional space time. This gives an obstacle to 
identifying the Feynman diagrams with smooth worldsheets of strings, which means that 
the bosonic model cannot be interpreted as a string theory. We expect that the situation 
changes drastically in the supersymmetric case. 

We defined a quantity representing the uncertainty of the space-time coordinates and 
showed that it is of the same order as the extent of space time, which means that a classical 
space-time picture is maximally broken. To all orders of the 1/D expansion, we found that 
correlation functions are given by an expansion in terms of 1/N 2 . This makes, e.g., N = 16 
sufficiently large to probe the large N behavior of the model even for large D. We pointed 
out that this is related to the maximal breakdown of the classical space-time picture. 

The absence of an SSB of the Lorentz invariance was shown both by the Monte Carlo 
simulation and by the 1/D expansion. 

All the above issues we addressed for the bosonic model should be addressed for the 
IIB matrix model as well. The results must reveal all the dynamical properties of the 
model related to the space-time structure and the worldsheet picture. We expect that the 
supersymmetry plays an essential role. We hope our findings for the bosonic model provide 
a helpful comparison when we investigate the supersymmetric case including the IIB matrix 
model. 
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Appendix: The algorithm for the Monte Carlo simula- 
tion 



In this appendix, we explain the algorithm we use for the Monte Carlo simulation of the 
bosonic model with the action 

S=-^$>(LVAf), (A.l) 

where j3 = jp. Since // ^ is, the action is quadratic with respect to each component, which 
means that we can update each component by generating gaussian variables using the heat 
bath algorithm. 

We first note that the trace can be written as 

tr(L4 M , A u ] 2 ) = 2{A^ ij {A v ) jk {A tl ) kl {A v ) li - 2(A fJ ,) ij (A u 2 ) jk (A IJ ,) ki . (A.2) 

Suppose we want to update A^. In the first term, all the A^ 2 components are coupled, 
while in the second term, only the components which include at least one common index are 
coupled. 

In order to simplify the algorithm, we employ the following trick. We first rewrite the 
action as 

S = -^E tr «A^} 2 ) + 2/3^tr(A M 2 ^ 2 ). (A.3) 

Here we introduce the auxiliary field Q^ v (1 < fi < v < D), which is N x N hermitian 
matrices, with the following action 

S' = ^X>W, " 2Q^G IW ) + 2{3 £tr(4X 2 )> ( A - 4 ) 
where is an hermitian matrix defined by 

G fiV = {A ll ,A v }. (A.5) 

Integrating out the auxiliary field Q^u, one reproduces the original action ( |A.3|) . 

Note that the simultaneous update of each of (Q^ij can be done easily by generating 
Gaussian variables. As for the update of A^, note that only the components which include 
at least one common index are coupled due to the second term of ( A.4 ). The A" diagonal 
components (A^a (i = 1, • • • , N) for each /x can be updated simultaneously. The y off- 
diagonal components (A M ) nina , (y4 M )„ 3n4 , (A^) nN _ iriN , where n 1 ,n 2 ,---,n N are different 
indices, can be updated simultaneously for each \x. 



39 



References 



[i] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7 



[9 

[10 

[11 
[12 
[13 
[14 



T. Eguchi and H. Kawai, Phys. Rev. Lett. 48 (1982) 1063. 

G. Bhanot, U. Heller and H. Neuberger, Phys. Lett. 113B (1982) 47. 

G. Parisi, Phys. Lett. 112B (1982) 463. 

D. Gross and Y. Kitazawa, Nucl. Phys. B206 (1982) 440. 

S. Das and S. Wadia, Phys. Lett. 117B (1982) 228. 

T. Banks, W. Fischler, S.H. Shenker and L. Susskind, Phys. Rev. D55 (1997) 5112, 
jhep-th/9610"04l . 



N. Ishibashi, H. Kawai, Y. Kitazawa and A. Tsuchiya, Nucl. Phys. B498 (1997) 467, 
|hep-th/9612~TT5 . 



M. Fukuma, H. Kawai, Y. Kitazawa and A. Tsuchiya, Nucl. Phys. B510 (1998) 158, 
liep-th/970512S . 



H. Aoki, S. Iso, H. Kawai, Y. Kitazawa and T. Tada, Prog. Theor. Phys. 99 (1998) 713, 
|hep-th/ 9802085 . 



W. Krauth and M. Staudacher, Phys. Lett. B435 (1998) 350, ^ep-th/9804TT9 



A. Gonzalez-Arroyo and M. Okawa, Phys. Rev. D27 (1983) 2397 

M. Okawa, Phys. Rev. Lett. D49 (1982) 705. 

J. Nishimura, Mod. Phys. Lett. All (1996) 3049, jhep-lat/9608TT9 



T. Yoneya, p. 419 in "Wandering in the Fields", eds. K. Kawarabayashi and A. Ukawa 
(World Scientific, 1987); p. 23 in "Quantum String Theory", eds. N. Kawamoto and T. 
Kugo (Springer, 1988); Mod. Phys. Lett. A4 (1989) 1587. 



[15] T. Nakajima and J. Nishimura, Nucl. Phys. B528 (1998) 355, Hiep-th/9802082 



40 



